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In the view of the quantum dynamic theory of measurement developed from the Hepp-Colemen 
(HC) model (K. Hepp, Hev.Phys.Acta, 45, 237 (1972)), the quantum decoherence in neutrino flavor 
oscillation caused by an environment surrounding neutrinos is generally considered in this paper. 
The Ellis, Hagelin, Nanopoulos and Srednicki (EHNS) mechanism for solving the solar neutrino 
problem can be comprehended in a framework of the ordinary quantum mechanics. In the weak- 
coupling limit, a microscopic model is proposed to describe the transition of two neutrino system 
from a pure state to a mixed state. It gives the modified formula of survival probability of neutrino 
oscillation with two additional time-dependent parameters. For specified environments, this result 
shows that the oscillating phenomena of neutrino still exist even without a mass difference in free 
neutrino. 
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It is well-known that environment may decohere a 
quantum system immersed in it, making this system to 
lose its coherence and to transit from a pure state to a 
mixed state [ 1]. Based on the quantum dynamic ap¬ 
proach of this environment influences [2], which was de¬ 
veloped from the Hepp-Colemen (HC) theory of quantum 
measurement and its clarified physical representation [3], 
the decoherence phenomenon is investigated for the neu¬ 
trino flavor oscillation [4] in this paper. 

Neutrino oscillation is a subject of much intense the¬ 
oretical and experimental research. In comparison with 
the prediction from the standard solar model(SSM), the 
deficit of solar neutrinos observed in the earth can prob¬ 
ably be explained by neutrino oscillations with a (mass)^ 
difference Ato^ ~ 10“®eV^ [4]. Two possible mechanisms 
were further presented to enhance the vacuum neutrino 
oscillations. As a matter enhancement effect of oscilla¬ 
tions, the Mikheyev-Smirnov-Wolfenstein (MSW) mech¬ 
anism [ 5], was invoked when the neutrinos pass through 
the sun adiabatically and non-adiabatically [6]. Initialed 
by Ellis, Hagelin, Nanopoulos and Srednick (EHNS) in 
1983 [7], there was another mechanism ( now called 
EHNS mechanism) to modify neutrino oscillations. Sim¬ 
ilar to the quantum decoherence in generic problem of 
environment influence [1-3], the EHNS mechanism makes 
a transition of neutrino from a pure state to a mixed state 
due to some unknown couplings of neutrinos to the envi¬ 
ronment surrounding them. Such an environment may be 
the background field caused by Hawking s quantum grav¬ 
ity effect, like the evaporation of black holes [8]. With 
irreversible elements characterized by phenomenological 
parameters, the original EHNS mechanism suggested an 
effective von Neumann equation of density matrix. Its 
solutions modify the survival probability of neutrino os¬ 
cillation . Following this work some significant efforts 
were made by different authors in preserving the linear¬ 
ity, locality, conservation of probability and unitarity [9]. 
Notice that the EHNS parameters were determined to 
some extent either by analyzing the experimental data of 
CP violation in system [10] or based on further 


theoretical considerations, such as the string theory with 
an Einstein-Yang-Mills black hole in four dimensions [11]. 

These investigations reflect the successful aspects of 
the EHNS mechanism in solving the solar neutrino prob¬ 
lem. However, as the phenomenological motion equa¬ 
tion of density matrix concerns quantum irreversible pro¬ 
cess, it apparently violates quantum mechanics since the 
Schroedinger equation is reversible with time-reversal 
symmetry. As a correct physical theory to solve the solar 
neutrino problem, the EHNS mechanism has to be consis¬ 
tent with quantum mechanics. To remove this confliction 
in this paper we show that rationality of the EHNS mech¬ 
anism is related to the reduction of quantum mechanics 
for the “universe” , a total system formed by the neu¬ 
trino system plus environment. The physical essence of 
our observation is that the couplings of neutrino to the 
EHNS mechanism results in quantum decoherence of the 
flavor eigenstates of neutrino, which are coherent super¬ 
positions of neutrino mass eigenstates if there is not an 
external influence. Mathematically, this physical process 
is described by a reduced density matrixes of neutrino 
system by tracing out over the variables of the environ¬ 
ment. The calculation of the reduced density matrixes 
explicitly leads to a modified formula of survival proba¬ 
bility of neutrino oscillation. It not only introduces ex¬ 
tra parameters to neutrino oscillation, but also implies 
a novel dynamic effect that the oscillating phenomena of 
neutrino may exists even without a mass difference in 
free neutrino. 

Methodologically, the starting point of the present 
study is the quantum dynamic approach [2] for decoher¬ 
ence problem(or wave function collapse, WFC) in quan¬ 
tum measurement. We recall that, in the traditional the¬ 
ory of quantum measurement , the WFC postulate is 
only an extra postulate added to the ordinary quantum 
mechanics. Under this postulate, once we measure an ob¬ 
servable and obtain a definite value ak the state of system 
must collapse into the corresponding eigen-state \k) from 
a coherent superposition \(j)) = '^}^Ck\k){k\. Through 
density matrix this process is described by a quantum 
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decoherence process p = |^)((/'| p = 
from a pure state to a mixed state. In quantum dynamic 
approach, both the measured system and the measuring 
apparatus obey Schroedinger equation and the dynamic 
evolution governed by their interaction is supposed to 
result in quantum decoherence for the reduced dynamics 
of the measured system under certain conditions, such as 
the macroscopic limit that the measuring apparatus con¬ 
tains an infinite number of subsystems. Notice that the 
environment coupling to neutrino usually can be consid¬ 
ered as a bath of infinite degrees of freedom. For quan¬ 
tum decoherence, there is a strong resemblance between 
the neutrino system plus the environment and the mea¬ 
sured system plus the measuring apparatus. So we can 
stndy quantum decoherence effect in neutrino oscillation 
by making use of the quantum dynamic approach devel¬ 
oped in treating quantum measurement. In the following 
our discussion will focus on the case of two generations, 
but it can be easily generalized to the case of three gen¬ 
erations. 

The neutrino weak interaction (flavor) eigenstates Ve 
and do not coincide with their mass eigenstates i>i and 
1^2 with masses mi and m 2 . Through a vacuum mixing 
angle 6, the neutrinos mixture is described as a two di¬ 
mensional rotation 

t'e = cos 0-I-sin0 (I) 

Vfj, = —vi sin 9 + 1/2 cos 9 

Without the interaction between neutrinos and other ob¬ 
jects, the vacuum neutrino oscillation result from free 
Hamiltonian 

Hn = Eii/ii/i + E2V2V2 ( 2 ) 

As a result, the survival probability that an electron- 
neutrino remain as itself is: 

(t) cs 1 - sin^ {29) sin^ (3) 

Here 5E = E 2 — Ei, Sm^ = m^ — mf and we have consid¬ 
ered masses of neutrinos are so small that their energy 
and momentum are very close. The above well-known 
result was obtained completely from the coherent super¬ 
position Eq.(I) of pure states 1/1 and 1/2 and it will be 
certainly modified if the mixed states of neutrinos were 
introduced from fluctuating interactions. 

Now, we consider the coupling of neutrinos to environ¬ 
ment, which may be a background provided by the cer¬ 
tain gravity effect such as the Hawking’s evaporation of 
black holes. In general, without referring to the concrete 
forms of the environment Hamiltonian iJc (q) of dynamic 
variables q , we can roughly write down the interaction 
Hamiltonian in the weakly -coupling limit 

Hj = [\iDiVi + \2V2V2]q (4) 


which is linear for both the density projections Pii/i{i = 
1,2) of neutrinos and the environment variable q. The 
coupling constants Ai {i = 1 . 2 )satisfying Ai 7 ^ A 2 im¬ 
plies that the couplings to environment have different 
strengths for different neutrinos i/i and 1 / 2 - For instance, 
in the quantum gravity background, perhaps there is 
Ai oc mf. Up to the first order perturbation for quan¬ 
tum dissipation problem , Caldeira and Leggett [12] have 
pointed out that, any environment weakly coupling to a 
system may be approximated by a bath of oscillators with 
linear interaction in the weakly coupling condition that 
“each environmental degree of freedom is only weakly 
perturbed by its interaction with the system”. For both 
quantum decoherence [2,14]as well as quantum dissipa¬ 
tion [12,13] , we have modeled a general environment 
with the linear coupling Hamiltonians similar to Eq.(4) 
and discussed the universality of such modeling scheme. 

Because the total system, governed by the total Hamil¬ 
tonian iJ = Hn + Hq -f Ell 1 is closed, the quantum me¬ 
chanics can work well to describe its dynamic evolution. 
Suppose the neutrino is initially in a flavor state Ve with 
a pure density matrix 

cos^ 9 cos 0 sin 0 A , , 

cos 9 sin 9 sin^ ^ / 

and the environment in a general initial state of den¬ 
sity matrix pc,the initial condition for the total system 
should be : 

p{t = 0) = (g) pG (6) 

By considering the quantum non-demolition features of 
the interaction , i.e., [Hi, iL„] = 0, [Hi, Hq] ^ 0,the total 
density matrix can be directly computed as 

p{t) = cos^ 9 ■ vivi (g) Ui{t)pGUl{t)+ 

+ SITL^ 9 ■ V2V2®U2{t)pGUl{t)] (7) 

-I-cos0sin0 • V 1 V 2 (g) Ui{t)pGU\{t) ey+){i6Et) + H.c. 
where 

Ui{t) = exp{-i[HG{q) + Ai^jt} (i = 1,2) ( 8 ) 

for i = 1 , 2 are the effective evolution operators of envi¬ 
ronment. They show that the back-actions on environ¬ 
ment exerted by 1/1 and 1/2 respectively are different for 
different coupling strengths \i{i = 1,2). Namely, the en¬ 
vironment can distinguishes different mass eigenstates of 
neutrino. 

As our main interest is only in the neutrino oscillation 
rather than the motion of the environment, we trace out 
over the variable of environment to obtain the reduced 
density matrix of the neutrino system. It is 

Pn{t) = cos^ 9 ■ Dii/i + + sin^ 9 ■ V 2 V 2 
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+ F{t) cos 0 sin0 • i'ii '2 exp(i ) + H.c. (9) 

2E 

of the neutrino system. Here, the decoherence factor 

Fit)=Trg[U^{t)U^{t)pG] ( 10 ) 

characterizes the extent that the nentrino system loses its 
coherence. The reduced density matrix Pn(t) determines 
the survival probability that a electron-neutrino remain 
in itself is 

= Trn[UeVePn{t)] 

= 1 — i sin^(20){l — cos[6Et + a{t)]\F{t)\} (11) 

where the real number a{t) is the phase of F{t). It is 
easy to prove that |F(t)| is less than 1 since Tr[pG] = 1 
and Ui{t) (i=l,2) are unitary. In the above concise for¬ 
mula ( 11 ) for neutrino oscillation, the influences of envi¬ 
ronment are simply summed up to a decoherence factor 
T'(t), which is time-dependent while energy-independent. 
If the difference Ai — A 2 in couplings of vi and to envi¬ 
ronment does not depend on (5m, the difference in their 
masses, the result ( 11 ) directly shows a novel fact that, 
even neutrino is massless or electron- and muon- neutri¬ 
nos have a same mass, it is possible to see the oscilla- 
tional behavior between electron- and muon- neutrino. 
Here, the neutrino oscillation is caused by the phase a{t) 
of the decoherence factor, that is 

= 1 - ^sin2(26»)(l -cos[a(t)]|F(t)|) ( 12 ) 

As the above derived temporal phenomenon is indepen¬ 
dent of the energy of neutrino flux, it can not solve the 
present solar neutrino problem with energy -dependence 
spectrum . 

To get one step further we consider two extreme cases. 
When |T'(t)| = 1, or F{t) = exp[iQ;(f)] for real a, the 
system still in a pure state 

Pn{t) = AAAt) : 


(y. Ck,(t] 

At) = vi exp(iifit — i—) cosO + V 2 exp(iif 2 t + 

(13) 


In this case the influences of environment on neutrino 
oscillation only result in a time-dependent phase shift 


= 1 - sin^(26')sin2[^^ -h (14) 

As a signature of the environment influence, this phase 
shift only affects the spectrum distribution P{E) = 
p^^^ue (t) in Ihs higher energy region. For an extreme 


case with E{t) = 0, the coherence in neutrino system 
is completely lost and the neutrinos transit into a com¬ 
pletely mixed state 

Pn{t) = cos^ 6 ■ Dii/i -I- -I- sin^ 9 ■ V2i’2 (15) 

with vanishing off-diagonal elements. It leads to a con¬ 
stant transition probability^ sin^( 20 ) from Vf. to and 
the corresponding survival probability that Vg remain in 
itself is 1 — 1 sin^(20). Notice that these results only de¬ 
pend on the vacuum mixing angle 9 and do not fit the 
neutrino spectrum varying with the energy of neutrino. 

To get the idea of explicit form of F{t) we illustrate the 
determination of decoherence factor in the harmonic os¬ 
cillator model of environment [3,14]. Let and aj be the 
creation and annihilation operators for the j’th harmonic 
oscillator in environment. The Hamiltonian of environ¬ 
ment takes the form Ft = X(j=i TtuJjO^a and its interac¬ 
tion with neutrino can be modeled as a linear coupling: 

N 

Hi = [Aii^i^i -I- X2p2i'2\ ^ ^9jiaf + aj) (16) 

i=i 

Let the environment be initially in the vacuum state 
| 0 )e = jOi)® | 02 )( 8 )---( 8 ) [Oat) where | 0 j)is the ground state 
of the j’th single harmonic oscillator, the corresponding 
decoherence factor 


N N 

F{t) = F{N,t) = n {0j\U2\t)U({t)\0A = n m) 

i=i i=i 

can be obtained by solving the Schroedinger equations of 
UAA {i = 1)2) governed by the Hamiltonians of forced 
harmonic oscillators Hf = TiuijO^aj +Xigj{a'^ +aj). The 
result is 


|F(t)|=exp[- 


N 

E 

i=i 


2gASXr^ 


■ sm 




LOG 


a{t) = 




sin(u;jt) 






UJi 


(17) 


where 5A = A 2 — Ai,(5A^ = A| — A^.The decoherence 
time is decided by the norm part of F{N,t), which is 
the same as that from the two level subsystems model 
of environment in the weakly coupling limit [2]. Espe¬ 
cially, when the environment consists of identical parti¬ 
cles, we have loj = to and pj = g. In the van Hove limit of 
weakly-coupling , i.e., g ^ 0,N ^ 00 , but gVN ^a 
finite constant G,a simple result for decoherence fac¬ 
tor is obtained as |E(t)| = exp[—((5A)^^^ sin^(^)] and 
a{t) = 6X'^[A^t + £2llEfl]_ Then, an explicit formula of 
environment-modifying neutrino oscillation is obtained 
as 
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, G^6X^. 

{l-cos[(—+ ^^)i- 


4:E 


2w2 


(5A 


2sin{ujt) 


2u! 


exp[-((5A)2^sin2(y)]} (18) 

In the low-frequency case of environment (lu —>■ 0),denote 
P = A(5A)^G^, we have an exponential decay |F(t)| = 
exp[—/3t] and a{t) = So 


. (t ) ^ 1 — - sin^ (20) X 


{l-cos[(^-f y)t]} xexp[-/3t]. (19) 

In conclusion, we have considered the quantum deco¬ 
herence problem in neutrino flavor oscillation based on 
our dynamic approach for quantum measurement. The 
phenomenon of transition of neutrinos from a pure state 
to a mixed state, similar to EHNS mechanism, has been 
understood to some extent in the view of the ordinary 
quantum mechanics. Our study leads to the modified 
neutrino oscillations with two additional time-dependent 
parameters. For specified environments, they show that 
the oscillating phenomena of neutrino still exist even 
without a mass difference in free neutrinos. However, 
we have not considered yet the influences of quantum 
decoherence on the MSW mechanism and the quantum 
dissipation effects due to the couplings to environment 
in off-diagonal form [gviV 2 -\- H.c\q [12, 13]. These effects 
perhaps are crucial to the final solution to solar neutrino 
problem, and the method used by this paper can be gen¬ 
eralized to handle them without any difficult in princi¬ 
ple. The main problem in present is that the explicit 
calculation of the decoherence factor F{t) depends the 
details of the environment and its couplings to neutrinos 
in certain extent, but up to date we have not a complete 
knowledge yet about a complex environment. If the en¬ 
vironment is provided by the Hawking’s evaporation of 
black holes in gravity field, there is not yet an universally 
acknowledged scheme to quantize gravity. Regardless of 
this problem, the present paper at least suggests a general 
rule to microscopically analyze the decoherence modifi¬ 
cation of neutrino oscillation. Once the Hamiltonians of 
the environment and its coupling to neutrinos, even in 
the weakly-coupling case, were provided the decoherence 
factor can be explicitly calculated to fit the datum in the 
solar neutrino problem. 
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